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MODIFICATION OF AN IMPULSE -FACTORING ORBITAL-TRANSFER 
TECHNIQUE TO ACCOUNT FOR ORBIT-DETERMINATION 
AND MANEUVER-EXECUTION ERRORS 

By James F. Kibler, Richard N. Green, 

George R. Young, and M. G. Kelly 
Langley Research Center 

SUMMARY 

A method has previously been developed to satisfy terminal rendezvous and inter- 
mediate timing constraints for planetary missions involving orbital operations. The 
method uses impulse factoring in which a two-impulse transfer is divided into three or 
four impulses which add one or two intermediate orbits. The periods of the intermediate 
orbits and the number of revolutions in each orbit are varied to satisfy timing constraints. 
In this paper, techniques are developed to retarget the orbital transfer in the presence of 
orbit-determination and maneuver-execution errors. Sample results indicate that the 
nominal transfer can be retargeted with little change in either the magnitude (AV) or 
location of the individual impulses. Additionally, the total AV required for the retar- 
geted transfer is little different from that required for the nominal transfer. A digital 
computer program developed to implement the techniques is described in the appendix. 

INTRODUCTION 

The objective of one type of planetary mission might be to land a package of instru- 
ments at a preselected site on the surface of a planet to study its physical characteristics. 
The Viking-Mars program is a mission designed to achieve such an objective, and it will 
be used as an example throughout this paper. To accomplish the goal of surface recon- 
naissance, the Viking mission plan is to insert a lander -orbit er combination into orbit 
about Mars. The orbiter engine then makes up to four maneuvers (assumed here to be 
impulsive) to correct errors made in the Mars-orbit-insertion (MOI) maneuver, position 
the vehicle to allow reconnaissance of the landing site, and synchronize the spacecraft 
over the landing site in the proper position to make the deorbit maneuver. The lander 
then separates from the orbiter, deorbits, and lands while the orbiter maintains a data 
relay link to Earth. The sequence of orbital maneuvers may be characterized as a time- 
fixed rendezvous from a specified initial orbit to a final orbit which has a synchronous 

period P , a periapsis radius r which does not violate planetary quarantine 
s p 



constraints, and a particular true anomaly fpgj^ directly over the landing site. The 
value of fp ER is defined by a study of the deorbit sequence. The landing-site recon- 
naissance may be interpreted as an intermediate timing constraint on the rendezvous 
sequence. 


One method of determining the proper sequence of orbital maneuvers is reported 
in reference 1. In this method, the problem is separated into two parts. First, a two- 
impulse maneuver is found which transfers from the post-MOI orbit (initial elliptical 
orbit) to a final orbit which satisfies the geometry constraints ^P g , Tp, and fp ER over 
the landing-site latitude 6p E p)* Then, one or both of the two impulses are factored 
into two or three collinear parts. For example, for the first impulse, 


AVj = qAVj + (1 - q)AV^ 

By applying the first part of the impulse (qAvJ, waiting one or more revolutions, and then 
applying the remainder of the impulse ((1 - q)AVjJ at the same point, there is no net 
change in the resulting geometry. However, qAV^ places the spacecraft in an inter- 
mediate orbit whose period is generally different from those of the two original orbits. 
This period difference and the number of revolutions in each orbit may be varied to 
satisfy the timing constraints. If the total number of impulses is limited to four, five 
different maneuver sequences arise. Two are three-impulse solutions (factor either the 
first or second impulse), and three are four-impulse solutions (factor each impulse, or 
factor the first or second impulse twice). 

The orbital maneuver sequences determined by the impulse-factoring method do not 
consider the effects of errors in the determination of the orbital elements nor of errors 
in the execution of each of the maneuvers. Due to errors in the determination of a space- 
craft trajectory from observation data, only an estimate of the actual trajectory is avail- 
able. A maneuver can be computed based on the estimate of the trajectory, but in reality 
it is applied to the actual trajectory. Further errors are introduced when the maneuver 
is applied since the maneuver cannot be executed precisely. One way to study the effects 
of these error sources is by a Monte Carlo analysis. In this method, the distribution of 
each error source is assumed. A random sample is obtained from each error distribu- 
tion and the mission is simulated. By repeatedly simulating the mission with different 
random samples, statistics may be accumulated which represent the random process. 

The procedures described in this paper represent a single simulation in such a Monte 
Carlo analysis. 

Only four-impulse transfers from the post-MOI orbit to the final orbit are consid- 
ered. Since the addition of errors at any point in the sequence of maneuvers invalidates 
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the remainder of the original no -error sequence, the remaining maneuvers must be 
retargeted from the new orbit to the given final orbit. The rationale for the three-, two- 
and one-impulse retargeters is presented in this paper. Sample results are given for 
retargeting a typical Viking-type mission. The appendix contains a brief description of 
the digital computer program which has been developed to implement the impulse fac- 
toring and retargeting technique. 


SYMBOLS 


f 

f PER 

F 

G 

i 

I,J,K,L,M 


k 

_P_ 


P 


s 


q,r 


*A 

t miss 


true anomaly, deg 

true anomaly of a point in orbit directly over the landing site, deg 
cost function 

timing error in final orbit, t^ - t D , hr 
index defining reconnaissance orbit 

integral revolutions in initial, first factored, transfer, second factored, and 
final orbits, respectively ' 

weighting factor 

.or.bital.period, hr _ .... 

synchronous orbital period, hr 
velocity factors 
periapsis radius, km 

actual time of arrival over landing site, hr 
desired time of arrival over landing site, hr 
target miss time, hr 
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Atf time from entrance to final orbit to fp£R> hr 

AV vector describing an impulsive velocity change, m/s 

AV magnitude of AV, m/ s 

AVg AV required for the two-impulse transfer which satisfies the geometry con- 

straints on the final orbit AV^ + AVg, m/s 

6 per landing-site latitude, deg 

Subscripts: 


1,2 refer to the first and second impulses of the geometry solution 

f final orbit 

Superscripts: 

r,h,n radial, horizontal, and normal components of AV, respectively 

ANALYSIS 


Nominal Orbital Transfer 

After the initial Mars orbit is established, an impulse-factoring technique can be 
used to find a sequence of four impulsive maneuvers which satisfy the following con- 
straints on the Viking mission: 

1. Transfer from the known initial orbit to a final orbit which has the correct 
period, periapsis radius, and true anomaly over the landing-site latitude. 

2. Allow reconnaissance of the landing site on a particular revolution after 
insertion. 

3. Synchronize the spacecraft with the landing site to allow the deorbit maneuver. 

The first step of the impulse-factoring technique requires obtaining a "geometry 
solution" which satisfies the first constraint. This solution consists of two impulses: 
the first, AVp places the spacecraft into a transfer orbit; the second, AVg, places the 
spacecraft into the geometrically correct final orbit. 
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In order to satisfy the second and third constraints, the two impulses may be fac- 
tored such that 


AV t = qAV 1 + (1 - q)AV 1 
AV 2 = rAV 2 + (1 - r)AV 2 


( 1 ) 


After I revolutions in the initial orbit the first "factor" of AVp qAV^, is applied. 
This impulse places the spacecraft into an intermediate orbit between the initial and 
transfer orbits. Then, after J complete revolutions in this intermediate orbit, the 
spacecraft is placed into the transfer orbit by applying the second factor, (1 - q)AVj, 
at the same point. Similarly, after K revolutions in the transfer orbit rAV 2 is 
applied, placing the spacecraft into a second intermediate orbit. After L complete 
revolutions in the second intermediate orbit the last impulse, (1 - r)AV 2 , places the 
spacecraft into the final orbit. 

Since the intermediate orbits have periods which in general are different from 
those of the initial, transfer, and final orbits, the time spent in the two intermediate 
orbits serves to alter the timing of the geometry solution. By selecting proper values 
for q, r, I, J, K, and L, the reconnaissance and rendezvous constraints can be 
satisfied. 

Two other four-impulse schemes are 


AVj^ = qAV^ + r(l - q)AV 1 + (1 - r)(l - q)AV 1 


av 2 = av 2 


where AVj is factored twice and AV 2 is applied in full; and 


AVj = AV t 


AV 2 = qAV 2 + r(l - q)AV 2 + (1 - r)(l - q)AV 2 


where AV^ is applied in full and AV 2 is factored twice. 
Two three-impulse schemes are 


AV t = qAV t + (1 - q)AV 1 l 
A V 2 = A V 2 } 


( 2 ) 


( 3 ) 


(4) 
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where AVj is factored once and AV 2 is applied in full; and 


AV 1 = AV t 


AV 2 = qAV 2 + (1 - q)AV 2 | 


( 5 ) 


where AV^ is applied in full and AV 2 is factored once. 

The technique for selecting an optimum set of values for q, r, I, J, K, and L 
for each of the factoring solutions is the subject of reference 2. There are many possi- 
bilities which must be investigated, but for the present paper it is assumed that an opti- 
mum set has been found. This set determines the nominal four-impulse transfer 
maneuver to which orbit-determination and maneuver -execution errors are applied. 

One way to display the sequence of maneuvers is illustrated in sketch (a). Target 
miss time is plotted versus revolutions from MOI. The target miss time is evaluated 



Sketch (a). - Nominal targeting sequence. 

once each orbit as the spacecraft crosses the landing-site latitude. It is the time 
required for the landing site to rotate to a point directly beneath the plane of the orbit. 
As an aid to visualizing the maneuver sequence, the target-miss -time evaluation points 
are connected by straight lines. As the periods of the orbits change with each impulse, 
the slope of the connecting line is changed. For the Viking mission, zero slope repre- 
sents a synchronous period while positive and negative slopes represent super- 
synchronous and sub -synchronous periods, respectively. The target miss time must be 
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zero both for the reconnaissance pass and when the spacecraft reaches the final orbit. 

The final orbit should have a synchronous period. By varying the factors and the number 
of revolutions in each orbit, the timing constraints can be satisfied. 

The Overall Retargeting Technique 

Although a nominal orbital transfer (here restricted to four impulses) has been 
established, the problem is not yet solved. Since the nominal four-impulse transfer can- 
not be performed perfectly because of orbit-determination and maneuver -execution 
errors, a method must be developed to retarget the individual impulses to correct the 
final target parameters within acceptable tolerances. The parameters which must be 
controlled on the final orbit for the Viking mission are periapsis altitude, period, target 
miss time, and fpgp^ over the landing-site latitude. A lower bound is placed on the 
periapsis altitude to satisfy the Mars quarantine constraint. The errors in the final- 
orbit period and target miss time are constrained since the spacecraft is required to be 
in a synchronous orbit directly over the landing site. The purpose of this requirement 
on the Viking mission is to allow the lander to deorbit to the landing site while the 
orbiter maintains a communications link with the Earth-based tracking stations. 

The overall technique for determining the effect of maneuver -execution and orbit- 
determination errors on the impulse sequence and final target parameters is as follows: 

1. Generate a nominal four-impulse transfer from the actual initial orbit to the 
final orbit using the impulse-factoring technique to satisfy all constraints. The result is 
a sequence of four impulses applied at particular points in each orbit. 

2. Perturb the actual initial orbit with orbit-determination errors, obtaining an esti- 
mate of the initial orbit. 

3. Retarget the predicted four-impulse transfer from the estimated initial orbit to 
the known final orbit in order to remove accumulated errors. 

4. Add execution errors to the first impulse of the retargeted four-impulse solution 
and apply this first impulse to the actual initial orbit. This step results in an actual sec- 
ond orbit. 

5. Add orbit-determination errors to the actual second orbit, obtaining an estimate 
of the second orbit. 

6. Using the estimated second orbit as an initial orbit, retarget a three-impulse 
transfer to the known final orbit in order to remove accumulated errors. 

7. Repeat steps 4, 5, and 6 by applying errors to the first maneuver and first orbit 
in each sequence, thus retargeting two- and one-impulse transfers. The imperfect orbital 
transfers result in errors in the final orbit which cannot be removed completely. By 
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keeping track of the actual location, magnitude, and direction of each of the four applied 
maneuvers, the errors in the final-target conditions can be calculated. If the errors are 
within acceptable tolerances for a wide range of sample cases, the retargeting strategy 
is successful. 

A typical retargeting sequence is illustrated in sketch (b). For clarity, target miss 
times versus revolutions from MOI are plotted as straight lines. The solid lines indicate 
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Sketch (b). - Typical retargeting sequence. 

the intended variation in target miss time as each impulse is retargeted. The dashed line 
traces the actual target miss time throughout the sequence. The difference between the 
dotted and solid lines represents the period error of each orbit due to maneuver-execution 
and orbit-determination errors. 

The following sections present one way of formulating each of the retargeting 
methods. 

Four -impulse retargeting . - The four -impulse retargeting is a simple modification 
of the original four-impulse problem. Orbit -determination errors are added to the actual 
initial orbit to yield an estimate of the orbital elements. A new four-impulse sequence 
is generated based on the estimated initial orbit. Then, the first impulse of the new 
sequence is perturbed with execution errors and applied to the actual initial orbit to yield 
an actual second orbit. 

Three-impulse retargeting . - The three-impulse retargeting is formulated as a 
modification to the three-impulse transfers (eqs. (4) and (5)). The initial orbit is taken 
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to be the actual second orbit that results from applying the first impulse of the retargeted 
four -impulse transfer with execution errors imposed. An estimate of this actual orbit 
is obtained by adding orbit-determination errors. A small timing error is allowed in the 
equations for the reconnaissance condition to account for accumulated errors at the time 
of landing-site reconnaissance. The AV required for the three-impulse retargeted 
transfer is minimized by the adaptive creeping algorithm discussed in reference 3. The 
velocity factor is adjusted to remove errors in the final target parameters. 

Two-impulse retargeting . - The two-impulse retargeting requires a formulation 
different from those previously used. The actual second orbit of the retargeted three- 
impulse transfer is used as the initial orbit for this portion of the trim strategy. An 
estimate of the initial orbit is found by adding orbit-determination errors. Since two 
impulses are required to transfer from this estimated initial orbit to the given final 
orbit, the impulse factoring technique can no longer be used. The only timing constraint 
which must be satisfied is the target miss time in the final orbit since the reconnaissance 
pass has been made on a previous orbit. A penalty -function approach is used to satisfy 
the final timing constraint. In this approach, the original cost function AV is aug- 

D 

mented by a weighted function of the equality constraint to be satisfied. The augmented 
function F is then successively minimized with increasing weights to drive the equality 
constraint to its desired value. That is, F is minimized where 

F = AV + kG 2 
g 

It is desired to drive G to within some acceptable tolerance. By choosing a sufficiently 
large k and minimizing F, the timing constraint G will be satisfied at the expense of 
AV . Thus, the retargeted two-impulse transfer satisfies all constraints on the final 

o _ __ _ _ ^ 

orbit. 

One-impulse retargeting . - For the one-impulse retargeting, the initial orbit is 
taken to be the second actual orbit of the retargeted two-impulse transfer. An estimate 
of this initial orbit is obtained by adding orbit-determination errors. Since only one 
impulse remains for the transfer to the desired final orbit, not all of the elements of the 
final orbit can be controlled simultaneously. Thus, a tradeoff must be made between the 
final timing error and the period error in the final orbit to keep both parameters within 
acceptable tolerances. The method used to target the final impulse is as follows: 

1. Given the estimated elements of the fourth orbit and the desired final orbit along 
with the number of revolutions in each, a final target miss time t m . may be computed. 

2. The miss time is nominally corrected to zero by targeting to an asynchronous 
final period: 
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d _ u *miss 
f ' s aU 
M+ «i 

s 

The new period of the final orbit has the effect of apportioning the timing error over M + 
revolutions. 

3. Given the new final period, a new semimajor axis of the final orbit is calculated. 

4. By holding the eccentricity of the final orbit and true anomaly into the final orbit 
constant, a one-impulse transfer can be computed to target to the new semimajor axis on 
the final orbit (essentially the second impulse of the two-impulse transfer described in 
ref. 1). Since the true anomaly is held constant, At^ will change and a new computed 
t m iss not necessarily be zero. If t migg is not within bounds, the entire proce- 
dure can be repeated from step 2 above using the new At^. Since there is no flexibility 
in the one-impulse transfer, AV cannot be minimized. 

Application of Orbit-Determination and Maneuver -Execution Errors 

Each of the retargeting methods discussed in the previous sections is applied to an 
initial orbit which has been corrupted by orbit -determination errors. The errors are 
obtained by sampling from one of two multivariate distributions which are represented by 
satellite-knowledge covariance matrices. The errors are in the classical orbital ele- 
ments (semimajor axis, eccentricity, inclination, argument of periapsis, right ascension 
of the ascending node, and time of periapsis passage). One covariance matrix is valid 
at apoapsis and the other is valid at periapsis. The estimated true anomaly for the next 
maneuver determines which matrix is to be used. That is, if the maneuver true anomaly 
is between -90° and 90° the periapsis covariance matrix is used; otherwise, the apoapsis 
matrix is used. These matrices are input to the program described in the appendix. 

After retargeting methods are applied, the estimated sequence of maneuvers is 
known. However, before the next impulse is applied maneuver -execution errors are 
added. The following errors (assumed independent) are considered: accelerometer 

bias, accelerometer calibration, ignition timing error, and errors in two pointing angles 
(right ascension and declination of the fixed thrust vector) for each of the four trim 
maneuvers. These errors have the effect of perturbing the elements of the succeeding 
orbit in the maneuver sequence. The standard deviation of each error source is input to 
the program described in the appendix. The orbit-determination and maneuver-execution 
errors used for the following sample results are representative of those expected for the 
Viking mission. 
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SAMPLE RESULTS 


The computer program developed to implement the retargeting techniques is 
described in the appendix. Figures 1 and 2 are a listing and flow chart of the main pro- 
gram. The subroutines required are also described in the appendix. A sample Viking 
case was computed to illustrate the retargeting techniques described herein. 

Input for this case is shown in figure 3 and described in the appendix. The actual 
initial orbit is taken from the output of a program described in reference 4. On the final 
orbit, fp ER is -10.6°, the target latitude is 12. 182°, and the final orbital period is 
synchronous (24.623 hr). Most of the remaining input is program related and does not 
change from case to case (see the appendix). The nominal four-impulse maneuver 
sequence from the actual initial orbit to the required final orbit is given in figure 4. In 
figure 5, the actual initial orbit has been perturbed by orbit-determination errors and the 
sequence has been retargeted to the required final orbit. In figure 6, the second orbit of 
figure 5 has been perturbed by orbit-determination and maneuver-execution errors and 
retargeted to the required final orbit. In figure 7, errors have been added to the second 
orbit of figure 6 and the maneuvers retargeted. In figure 8, errors have been added to 
the second orbit of figure 7. Here, however, a slightly asynchronous orbit of 
24.62355 hours is required in order to minimize the final timing error. The following 
table summarizes the retargeted maneuvers: 


Maneuver 

Impulses 

Timing 

error, 

hr 

First 

Second 

Third 

i 

Fourth 

deg 

AV, 

m/s 

i 

f, 

deg 

AV, 

m/s 

f, - 

deg 

AV,. 

m/s 

f, 

deg 

A V 

m/s 

Nominal 

four-impulse 

110.4 

35.3 

112.7 

4.4 

113.0 

7.7 

-152.4 

n.i 

-0.005 

Retargeted 

four-impulse 

- - 1 

109.6 

35.6 

111.9 

5.2 

112.2 

6.9 

-156.6 

10.7 

-0.006 

Retargeted 
three -impulse 



109.2 

3.9 

109.4 

8.0 

-147.7 

11.8 

-0.0005 

Retargeted 

two-impulse 





110.3 

7.9 

-148.0 

12.5 

0.0 

Retargeted 
one -impulse 







-148.0 

11.5 

0.004 
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The small changes in the retargeted controls for this sample case demonstrate the sta- 
bility of the nominal maneuver sequence in the presence of orbit-determination and 
maneuver-execution errors. The summary of actual end conditions in figure 8 shows 
that the error in the final orbital period is 7 seconds, the error in final periapsis radius 
is 0. 75 km, and the actual final timing error is -5. 3 seconds. All of these errors are 
well within acceptable tolerances. The actual total AV required for the trim is 
58. 5 m/s as opposed to a nominal total AV of 58.4 m/s, indicating that no appreciable 
penalty is encountered to correct orbit-determination and maneuver-execution errors for 
the sample case. Experience has indicated similar behavior for a wide range of maneu- 
ver sequences. 


CONCLUDING REMARKS 

A method has been presented to retarget an impulse -factoring orbital transfer to 
account for orbit-determination and maneuver-execution errors. Sample results indi- 
cate that the errors can be eliminated at small cost in velocity change by the retargeting 
techniques developed here. Additionally, the original solution for the sample case is 
quite stable in the presence of errors since the retargeted solutions are close to the nom- 
inal both in magnitude and location of the impulses. A digital computer program has 
been developed at the Langley Research Center to implement the impulse-factoring and 
retargeting techniques. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 4, 1974. 
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APPENDIX 


PROGRAM DESCRIPTION 
General 

A digital computer program has been developed to implement the impulse-factoring 
and impulse-retargeting technique described in the text. The program is coded in 
FORTRAN for the CDC 6600 computer. All of the subroutines required for MITOP 
(Multiple Impulse Transfer Optimization Program) are stored on the Langley Research 
Center data cell. The core storage required is about 70 K octal locations. Running time 
per case ranges from 1 to 10 sec of CPU time depending on the types of options exercised. 
Almost all data transfer is through two common arrays (C and IC) for floating point and 
integer variables. The contents of these arrays are described in tables I and II. Data 
associated with applying and accumulating errors for the retargeting cases are transferred 
through labeled COMMON TXERR described in table m. The following sections describe 
the main program, the primary-level subroutines, the input, and the output options. 

Main Program 

A listing of the main program for MITOP is presented in figure 1. A schematic 
flow chart of the main program logic is given in figure 2. The main program is simply 
an executive routine to handle input, output, data transfer, and scheduling between the 
primary-level subroutines. The modular construction of the MITOP-logic flow facilitates 
modifications to the program without disturbing the underlying subroutines. 

- - — - Primary -Level- Subroutines— • - 

The primary-level subroutines are those called by the main program. These 
special purpose subroutines in turn call the other subroutines listed on the data cell. 

The following is a brief description of each of the primary subroutines: 

INITIAL This subroutine is called to initialize inputs or various parameters which 

need to be calculated only once in a run. It is called with a single integer 
argument (I): 

1=0 Standard or "canned" inputs are stored (see tables I and II). 

1=1 Calculates parameters which remain constant for a single case. 

1=2 Calculates parameters which remain constant for a single 

Monte Carlo sample (this option is called from the other 
primary -level subroutines). 
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APPENDIX - Continued 


STATX 


RITE 


ERRORS 


BURN 4, 
BURN 2, 


This subroutine calculates the mean and standard deviation of 
parameters of interest which are input in a vector. The first 
call initializes various storage locations and subsequent calls 
accumulate statistics. It is called only in the Monte Carlo 
mode. 

This subroutine contains the various output options which are 
chosen by a single integer argument. The options are 
described in the output section. 

This subroutine adds errors due to maneuver execution and orbit 
determination to the initial orbit involved in the current retar- 
geting. It is called with a single integer argument which 
denotes the number of the maneuver being performed. The 
first call is with a zero to add orbit-determination errors to 
the initial orbit. This subroutine calls SAMPLE which gener- 
ates the errors in orbital elements by sampling from covariance 
matrices which are input. In addition, SAMPLE keeps track of 
the actual performed maneuvers and computes the errors in the 
actual final target parameters. 


BURN 3, ^ 
BURN !• J 


These subroutines generate the optimal four-, three-, two-, and 
one-impulse transfers, given an initial orbit and the final tar- 
get parameters. The BURN 4 subroutine is called twice. The 
first call is to generate a nominal four-impulse control history. 
The second call of BURN 4 (only in the retargeting option) is 
made after errors have been added to the initial orbit. A new 
four-impulse control history is calculated based on the per- 
turbed initial orbit. (The same technique is used for BURN 3, 
BURN 2, and BURN 1.) These subroutines call several other 
special purpose subroutines. MODEL is the function evalua- 
tion routine which calculates the two-impulse geometry trans- 
fer (GEOM), and depending upon the mode, satisfies the timing 
constraints (TRIM) on every iteration, or for only the geometry 
minimum transfer. GPOP is the function minimization routine 
which uses an algorithm due to Rosenbrock (ref. 3). GUESS 
discretizes the independent-variable space to provide a good 
first guess for GPOP. RANGE determines if the reconnaiss- 
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APPENDIX - Continued 


ance constraints are satisfied and MISS determines if the final 
timing constraint is satisfied. 

STATOUT This subroutine computes and outputs statistics on variables of 

interest. It is called only in the Monte Carlo mode. The out- 
put consists of mean and standard deviation of each parameter, 
the covariance matrix which indicates the degree of correlation 
between parameters, and histograms which indicate the distri- 
bution of selected parameters. 

Input 

The primary method of input to MITOP is by the NAMELIST identified as CASE. 
Selected variables from the C- and IC -arrays and from the TXERR COMMON are input 
by name. Any other required input may be made directly to the arrays. A sample input 
is shown in figure 3. The input variables are equivalent to locations in the C- and 
IC -arrays as described in tables I and II. If only a single case is run, set NMC = 1 
and input the initial orbit in AO. If a Monte Carlo case is run, the individual samples 
are input by a tape (identified as TAPE 10) which contains many post-MOI orbits from 
program VEAMCOP (ref. 4). The tape is read as follows: 

First record Total number of samples on tape (NOMC). 

Second record Five elements of the nominal initial orbit (ANOM). 

Third and subsequent records Five elements of the current sample initial orbit 

(AO), true anomaly of cutoff of the MOI maneuver 
(FBO), initial time bias (TBIAS), and the 
MOI AV (DVMOI). 

Multiple cases may be run with little effort. Only variables which change between cases 
must be input. 


Output 

There are several output options available in MITOP. The output is controlled by 
subroutine RITE which is called at various points in the program. The subroutine has a 
single integer argument which controls the type of print at each call. Four input loca- 
tions — IC(7 to 10) — are available to store the integer control key. Each location con- 
trols the output option chosen as follows: 
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APPENDIX — Continued 


IC(7) After nominal four-impulse targeting. 

IC(8) Prior to each retargeting. 

IC(9) After each retargeting. 

IC(10) After each Monte Carlo sample. 

Each output control key may take values of 0 to 10. The various values produce the fol- 
lowing output: 

0 No output. 

1 Comprehensive output (explained below). 

2 Minimal output consisting of one line (explained below). 

3 Output the contents of the C- and IC- arrays. 

4 Error summary at end of Monte Carlo sample (explained below). 

5 1 + 3. 

6 1 + 4, 

7 2 + 3. 

8 2 + 4. ■ 

9 1 + 3 + 4. 

10 Available for user-supplied output. 

Comprehensive output . - An example of the comprehensive output mode is given in 
figure 4. The labels are described as follows: 

ORBIT NUMBER 1 to 5 refer to initial orbit through final orbit. 

A semimajor axis. 
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APPENDIX - Continued 


E 

I 

W 

o 

F-IN 

F-OUT 

P 

HP 

HA 

PERIOD 
DELT 
PER LAT 

RADIUS -IN 

RADIUS-OUT 

R-OUT 

V-OUT 

DV 

TIME/ORBIT 


eccentricity. 

inclination. 

argument of periapsis. 

right ascension of ascending node. 

true anomaly of entrance to orbit. 

true anomaly of exit from orbit. 

semilatus rectum. 

periapsis altitude. 

apoapsis altitude. 

period. 

time from F-IN to F-OUT. 

latitude directly beneath fp E ^ (sub-PER 
point). 

magnitude of orbital radius at F-IN. 

magnitude of orbital radius at F-OUT. 

X-, Y-, and Z-components of orbital radius 
at F-OUT. 

components of orbital velocity at F-OUT. 

components of AV required to transfer 
between successive orbits. 

total time spent on each orbit. 
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APPENDIX — Continued 


CUM TIME OUT 

cumulative time from F-IN on the initial orbit 
to F-OUT on each orbit. 

REVS/ORBIT 

number of revolutions on each orbit. 

CUM REV TO FOUT 

cumulative revolutions to F-OUT on each orbit. 

DELTA V 

magnitude of each maneuver. 

MOI DV 

Mars -orbit-insertion AV. 

TOTAL DELTA V 

sum of all maneuvers. 

RECONNAISSANCE 

the reconnaissance pass occurs on the orbit 


(OCCURS ON ORBIT i) identified by orbit number i. 

The following reconnaissance-pass variables are labeled at three times (closest slant 
range approach to the target, spacecraft at fp ER , and the spacecraft directly over the 
target latitude): 


REV 
T. A. 

PHOTO TIME 

TARGET RA 

SUB S/C DEC and RA 

TIME TO IMPULSE, LAST and NEXT 

LOOK ANGLE 

RANGE 

CENT-ANG 


revolution from F-EN on the initial orbit, 
true anomaly. 

time from F-IN on the initial orbit. 

right ascension of target at PHOTO TIME. 

declination and right ascension of the spacecraft. 

time between the last impulse and PHOTO TIME 
and between PHOTO TIME and next impulse. 

angle at target between spacecraft and local 
vertical. 

slant range between spacecraft and target, 
central angle between spacecraft and target. 
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APPENDIX - Continued 


The following variables refer to the final orbital alinement. 

TARGET DEC and RA declination and right ascension of target when 

the spacecraft is at fp^R* 

SUB-PER DEC and RA declination and right ascension of sub-PER 

point. 


ANGULAR MISS 


central angle between target and spacecraft at 


f PER- 

MISS TIME AT PER timing error between target and spacecraft at 

f P ER' 


Minimal output . - The minimal output consists of a single line of numbers which 
are listed from left to right as follows: 


1 

2, 3, 4, 5 
6 

7,8 

9 

10 
11 

12,13,14 

15 

16 

17 

18 


Monte Carlo sample number. 

integer revolutions on each orbit. 

type of trim solution. 

q and r factors applied to AV. 

AV required to satisfy geometry constraints, 

— km/s. -- • — - - • - — — - - - 

AV required for total trim, km/s. 

true anomaly out of the initial orbit (F-OUT), deg. 

components of AV^, km/s. 

true anomaly into the final orbit (F-EN), deg. 

slant range at reconnaissance, km. 

final miss time, hr. 

MOI AV, km/s. 
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APPENDIX - Concluded 


Error summary at end of Monte Carlo sample . - This output option consists of the 
following variables which summarize the actual end conditions for each sample: 


FINAL ORBIT 


a, e, i, w, fi, F-IN, F-OUT on the final orbit. 


PDF 


period of final orbit, hr. 


DPDF 


difference between actual and desired PDF, s. 


HPF 

DHPF 

JDLAND 

TBIASF 

DVTRMA 


periapsis altitude of final orbit, km. 

difference between actual and desired HPF, km. 

actual Julian date when spacecraft reaches 
fpER in final orbit, days. 

timing error in the final orbit, s. 

actual AV required for trim maneuvers, 
km/s. 


In addition to the above print options, the following is output at the end of a Monte Carlo 
case: 


(1) A summary of each sample consisting of case number; AV^OT A ^TRIM’ 

AVtotal’ P er i°d> periapsis -altitude, and timing errors in the final 
orbit. 

(2) The means and covariance matrix of a V M qj5 av geOMETRY’ AV TRIM ; 

AVtqtal’ P er ^> periapsis -altitude, and timing errors. 

(3) Histograms of AV MOI ; A V TRIM ; AV TOTAL ; 311(1 P eriod > periapsis- 

altitude, and timing errors. 

(4) A listing of the summary statistics ordered on av totAL' 
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TABLE I. - DESCRIPTION OF ELEMENTS IN C -ARRAY 


Location 

Variable name 

Stored value 

Description 

1 

U 

42 828. 4 

Gravitational constant for Mars, 

km3/s2. 

2 

RS 

3393:4 

Radius of Mars, km. 

3 

DR 

2tt/360 

Degrees to radians conversion factor. 

4 

RD 

3 60/ 277 

Radians to degrees conversion factor. 

5 

RHOMAX 

2000.0 

Maximum acceptable slant range from 
the landing site to the spacecraft at 
reconnaissance, km. 

6 

SCONV 


No longer used. 

7 

CONV 


No longer used. 

8 

TLAT 

-30.0 

Target latitude, deg. 

9 

STLAT 


Sine of target latitude. 

10 

CTLAT 


Cosine of target latitude. 

11 to 15 

STEPLIM 


No longer used. 

16 to 20 

STEPPAR 


No longer used. 

21 to 25 

XIPLOW 

0, 0, 0.03, 
0, 0 

Lower bounds on independent vari- 
ables (F-OUT, AVj, AVp AVj, 
F-In) for discretization routine 
(GUESS). 

26 to 30 

XIPSTEP 

15.0, 0, 0.01, 
0, 15.0 

Step sizes for independent variables 
(above) for subroutine GUESS, 
km/s and deg. 


For the following, O refers to initial, T refers to transfer, 
and F refers to the final orbits, respectively. 

31, 61, 91 

AO, AT, AF 


Semimajor axis, km. 

32, 62, 92 

EO, ET, EF | 


Eccentricity. 

33, 63, 93 

XIO, XIT, XIF 


Inclination, deg. 
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TABLE I. - DESCRIPTION OF ELEMENTS IN C-ARRAY - Continued 


Location 

Variable name 

Stored value 

Description 

34, 64, 94 

WO, WT, WF 


Argument of periapsis, deg. 

35, 65, 95 

OO, OT, OF 


Right ascension of the ascending 
node, deg. ) 

36, 66, 96 

FOlN, FTIN, 

o, - 

True anomaly of entrance to the 


FFIN 


orbit, deg. 

37, 67, 97 

FOOUT, FTOUT, 
FFOUT 

-, -, -10.6 

True anomaly of exit from the 
orbit - FFOUT is defined to 

be *PER’ deg ‘ 

38, 68, 98 

SIO, SIT, SIF 


Sine of inclination. 

39, 69, 99 

CIO, CIT, CIF 


Cosine of inclination. 

40, 70, 100 

SWO, SWT, SWF 


Sine of little omega. 

41, 71, 101 

CWO, CWT, CWF 


Cosine of little omega. 

42, 72, 102 

SOO, SOT, SOF 


Sine of capital omega. 

43, 73, 103 

COO, COT, COF 


Cosine of capital omega. 

44, 74, 104 

SFOIN, SFTIN, 
SFFIN 


Sine of true anomaly of entrance. 

45, 75, 105 

CFOIN, CFTIN, 


Cosine of true anomaly of entrance. 


CFFIN 



"46", 76, 106 

sfoout;sftout, 

SFFOUT 


Sine of true anomaly of exit. ' 

47, 77, 107 

C FOOUT, C FTOUT, 
C FFOUT 


Cosine of true anomaly of exit. 

48, 78, 108 

SLRO, SLRT, SLRF 


Semilatus rectum, km. 

49, 79, 109 

PO, PT, PF 

-, -, 24.623 

Period, hr. 

50, 80, 110 

DELFO, DELFT, 


Change in true anomaly from 


DELFF 


entrance to exit, deg. 

51, 81, 111 

DELTO, DELTT, 


Time from entrance to exit, hr. 

' 

DELTF 


: 

52, 82, 112 

HAO, HAT, HAF 


Apoapsis altitude, km. 
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TABLE I. - DESCRIPTION OF ELEMENTS IN C -ARRAY - Continued 


Location 

Variable name 


Description 

53, 83, 113 

HPO, HPT, HPF 

-, -, 1500.0 

Periapsis altitude, km. 

54, 84, 114 

TBO, TBT, TBF 


Time bias, hr. 

55, 85, 115 

RAOP, RATP, 
RAFP 

• - 

Right ascension of the reconnais- 
sance point, deg. 

56, 86, 116 

FOP, FTP, FFP 


True anomaly of the reconnais- 
sance point, deg. 

57, 87, 117 

TFOIN, TFTIN, 
TFFIN 


Time of entrance to orbit referenced to 
periapsis, hr. 

58, 88, 118 

TFOOUT, TFTOUT, 
TFFOUT 


Time of exit from orbit referenced to 
periapsis, hr. 

59, 89, 119 

TFOP, TFTP, 
TFFP 


Time of reconnaissance point refer- 
enced to periapsis. hr. 

60, 90, 120 



Not used. 

For the following, Q refers to the first factored orbit (q-orbit) 
and R refers to the second factored orbit (r-orbit). 

121, 141 

AQ, AR 


Semimajor axis, km. 

122, 142 

EQ, ER 


Eccentricity. 

123, 143 

XIQ, XIR 


Inclination, deg. 

124, 144 

WQ, WR 


Argument of periapsis, deg. 

125, 145 

OQ, OR 


Right ascension of ascending node, deg. 

126, 146 

FQ, FR 


True anomaly of maneuver point, deg. 

127, 147 

SIQ, SIR 


Sine of inclination. 

128, 148 

CIQ, CIR 


Cosine of inclination. 

129, 149 

SWQ, SWR 


Sine of little omega. 

130, 150 

CW Q, CWR 


Cosine of little omega. 

131, 151 

SOQ, SOR 


Sine of capital omega. 

132, 152 

COQ, COR 


Cosine of capital omega. 
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TABLE I. - DESCRIPTION OF ELEMENTS IN C-ARRAY - Continued 


Location 

Variable name 

Stored value 

Description 

.133, 153 

SFQ, SFR 


Sine of maneuver true anomaly. 

134, 154 

CFQ, CFR 


Cosine of maneuver true anomaly. 

135, 155 

PQ, PR 


Period, hr. 

136, 156 

HAQ, HAR 


Apoapsis altitude, km. 

137, 157 

HPQ, HPR 


Periapsis altitude, km. 

138, 158 

QV, RV 


Factors associated with splitting 
AV. 

139, 159 

QP, RP 


Factors associated with splitting 
the change in period between 
orbits . 

140, 160 



Not used. 

161 

DV1 


AV^ Magnitude of the first 
maneuver of the two-impulse 
transfer which satisfies all . 
geometry constraints, km/s. 

162 

DV2 


AV 2 Magnitude of the second 
maneuver, km/s. 

163 

DVGEOM 


1 AVg Sum of AV^ +AV 2 , 

km/s. . . — . . - - 

164 

DVTRIM 


Total AV required to satisfy 
geometry and timing constraints, 
km/s. 

165, 166, 167 

DV1R, DV1H, 
DV1N 


Radial, horizontal, and normal 
components of AV^, km/s. 

168, 169, 170 

H2X, H2Y, 
H2Z 


X-, Y-, and Z-components of the 
angular momentum vector 
associated with the second orbit, 
km^/s. 
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TABLE I. - DESCRIPTION OF ELEMENTS IN C-ARRAY - Continued 


Location 

Variable name 

Stored value 

Description 

171, 172, 173 

PX, PY, PZ 


X-, Y-, and Z-components of a unit 
vector pointing toward PER on 
the final orbit. 

174 to 182 

RXYZ 


3x3 array containing the coeffi- 
cients for a coordinate transfor- 
mation from the PQW to the 
XYZ system. 

183 

RATO 


Right ascension of the landing site 
when the spacecraft is at peri- 
apsis on the initial orbit, deg. 

184 

ANGRAT 

360.0/24.623 

Angular rotation rate of Mars, 
deg/hr. 

185 

SH 

1.0/3600.0 

Seconds to hours conversion factor. 

186 

PMISS 

0.5 

Maximum acceptable miss time at. . 
reconnaissance, hr. 

187 

TMIN 

0.5 

Minimum time between reconnais- 
sance and the closest maneuver, 
hr. 

188 

TBIAS 


Time required for the landing site 
to rotate from its position when 
the spacecraft is at periapsis on 
the initial orbit to a point beneath 
the sub -PER right ascension, hr. 

189 

DELPTO 


Difference in period from the initial 
to the transfer orbit PT-PO, hr. 

190 

DELPFT 


Difference in period from the trans- 
fer to the final orbit PF-PT, hr. 

191 to 195 

STEPITR 

5.0, 0.003, 0.003, 
0.003, 5.0 

Initial step sizes in the independent 
variables (above) for the minimi- 
zation routine, deg, km/s. 
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TABLE I. - DESCRIPTION OF ELEMENTS IN C- ARRAY - Concluded 


Location 

Variable name 

Stored value 

Description 

196 

DVMOI . 


AV required for the Mars orbit 
insertion maneuver, km/s. 

197 

RHOMISS 


Slant range at the reconnaissance 
point, km. 

198 

RAMISS 


Difference in right ascension 
between the loading site and the 
spacecraft at reconnaissance, deg. 

199 



Not used. 

200 

TMISS 


Miss time in final orbit, hr. 




TABLE EL - DESCRIPTION OF ELEMENTS IN IC- ARRAY 


Location 

Variable name 

Stored value 

Description 

1 

NMC 

25 

Number of Monte Carlo samples. (If 
set to 1, only a single case is run. 
If >1, sample input is by tape. ) 

2 

IPHOTO 

8 

Revolution on which reconnaissance 
takes place. 

3 to 6 

IMEN, JMIN, 
KMIN, LMIN 

1, 2, 1, 2 

Minimum number of revolutions on the 
initial, first-factored, transfer, and 
second-factored orbits, respectively. 

7 to 10 


0, 0, 0, 4 . 

Write keys which may take values from 
0 to 9 (explained in output section). 

11 


1 

Retargeting control key: if 0 no retar- 
geting; if 1 retarget each impulse. 

12 to 15 

IBEST, JBEST, 
•KBEST, LBEST 


Distribution of revolutions in each orbit 
for the most favorable total trim. 

16 

KKGEOM 


Not used. 

17 

KEYBEST 


Best type of trim solution. 

18 

MODE 

1 

Cost function to be minimized: if 1 
minimize AV^; if 2 minimize 

av total = AV g + P enalt y due to 

timing; if 3 minimize AV + KG^. 

o 

19 

IMC 

1 

Number of the current Monte Carlo 
sample. 

20 

NREV 

19 

Periapsis passage on which transfer 
must be complete. 

21 to 25 

NSTEP 

24, 1, 2, 1, 24 

Number of steps in each independent 
variable for suoroutine GUESS. 

26 

ISUM 


Number of integer revolutions in trim 
solution. 

27 

ITR 

2500 

Number of iterations for subroutine 
GPOP. 
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TABLE H. - DESCRIPTION OF ELEMENTS IN IC -ARRAY - Concluded 


Location 

Variable name 

Stored value 

" — "H 

Description 

28 

ITRIM 

0 

Key to choose factoring method: if 
0 factor periods; if 1 factor AV. 

29 

1HALF 

0 

Key to choose l-g- minimum revolu- 
tions on each orbit: if 0 no constraint; 
if 1 at least 1^ revolutions on initial 
and transfer orbits. 

30 

| 

. 


Not used. 

31 to 36 

IKEY1 to IKEY6 

1, 1, 1, 1, 1, 1 

Keys to choose types of trim solutions: 
if 0 skip solution; if 1 test solution. 

37 to 49 



Not used. 

50 

MBEST 

2 

Number of integer revolutions on final 
orbit. 
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TABLE m. - DESCRIPTION OF VARIABLES IN TXERR 


Variable 

Dimensions 

Description 

TJDCA 

1 

Julian date of periapsis passage on nominal 
approach hyperbola. 

KNSAT 

1 

Satellite knowledge flag: 0, use covariance of 
errors; 1, perfect knowledge. 

VCAL 

1 

Calibrated value of velocity counting 
accelerometer. 

IOSIG 

1 

Special flag to indicate use of 1-sigma level 
satellite orbit determination and maneuver - 
execution errors: 0, off; 1, on. 

SATKN 

6x6x2 

6x6 satellite-knowledge covariance matrices 
valid at periapsis, SATKN(i,j,l), and apo- 
apsis, SATKN(i,j,2). These matrices repre-, 
sent errors in the classical orbital elements, 
a, e, i, co, ip* 

OTCONT 

11 

Standard deviations of trim control parameters 

and spacecraft parameters cr,. , , <L , 

°r € t IGN 

a l 6 1 a 2 5 2 a 3 6 3 a A 6 4 
6 is accelerometer bias; € is acceler- 

ometer calibration error; is ignition 

timing error; a and 6 are pointing angle 

errors for each of a possible four trims. 

TRACT 

7x5 

Array of values of classical orbital elements 
(a, e, i, (jo , S2, F-IN, F-OUT) for each of 
five possible actual ellipses (Post-MOI, Post- 
OT1,. ... POST-OT4)! 

PDA 

5 

Period of each of the five actual ellipses men- 
tioned above. 
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TABLE IE. - DESCRIPTION OF VARIABLES IN TXERR - Concluded 


Variable 

Dimensions 

Description 

TJDPNA 

6 

Julian date of initial periapsis passage for each of 
the five actual ellipses and the Julian date of 
periapsis passage on the orbit following the orbit 
maneuver. 

HPF 

1 

Periapsis altitude of the actual final orbit. 

DVTRMA 

1 

Total actual AV. 

DTAUF 

1 

Difference between actual final period and 
synchronous. 

DHPF 

1 

Difference between actual final periapsis altitude 
and 1500 km. 

TBIASF 

1 

Error in terms of Mars rotation time between the 
sublongitude of the PER point and the desired 
landing site. 
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PROGRAM MITOP ( I NPUT * OUTPUT * T APE 5= I NPUT ♦ T APE6=0UTPtJT . T APE 1 0 ) 

DIMENSION AO 15) .XIPLOW45) »XIPSTEP<5) .NSTEP ( 5 ) , XMEAN <7) .XCOV (7. 7) , 

• X I N 4 7 ) .ANOM (8) .XI P (5) .STEP (5) .SAVE <6. 100) . KE YOPT (6) 

S.CSAVE(200) . ICSAVE (50) 

COMMON C(200) .ICISO) 

COMMON/TXERR/TJDCA.KNSAT.VCAL. IOSIG .SATKN (6 ♦ 6. 2) .OTCONT I1I)« 

• TKACT (7.5) .PDA (5) .TJDPNA(6) »HPf .DVTRMA .DTAUF ,OHPF ♦ TBI ASF 

NAMEL I ST/CASE/AO. PER* TL AT. TB I AS. DVMOI.TMJN.NREV.IPHOTO.IMIN.JMIN. 
*KMIN*LMIN.NMC.KEVW1 .key W 2.C0NV .XIPLOto.XIPSTEP.NSTEP. MODE. I TP. IHALF 
*,KEYOPT.C.IC 

. .TJDCA.KNSAT. IOSIG. SATKN, OTCONT. VCAL 
EQUIVALENCE (A041) .C4311 ) . t PER. C 1 97)) . ITLAT.C(fi)) . 4TBI AS.C 4 188) ) . I 
®OVMOI .CC196) ) • (TMIN.CU 87) > • ( NREV , IC { 20 ) ) , (I PHOTO. IC (2) > » < IMIN. I C ( 
*3) ) . (JMIN.ICIA) > , (KMIN.ICI5) > ♦ (LMIN»IC(6) ) . (NMC.IC(l) ) . (KEYW1 . IC I 7 
®) ) . 4KEYW2.ICI8) ) . (CONV »C ( 7) ) . (XIPLOWU) ,C(21> ) . (XIPsTEP(l) ,C(26> ) . 
• < NSTEP ( 1 ) ,IC(2I)).(ITR.IC(27) ) , I IMC. IC 4 19) ) , (MODE. IC 1 18) ) . (FOOUT.C 

• (37) ) . (DV1R.C(165M . (0V1H.CU66) ) , (DV1N.C(167) ) . (FFIN.C(96) ) . (DVGE 
•OH. C (163) ) , (DVTRIM.CU6A) ) , ( I TR I M , IC (20 )) , < IHALF , I C ( 29) ) . (KEYOPTU 
«) . IC ( 31 » ) 

CALL INITIAL(O) 

MBEST-IC (50) =2 

1 CONTINUE 

READ (5. CASE) 1 IF(E0F,5> 2.3 

2 STOP 

3 CONTINUE 
I C (9) =2 
IC (9) =1 

WR I TE (6 .CASE ) 

CALL INITIALU) 

DO 160 1=1*200 

160 CSAVE(I)=C(I) 

DO 161 1=1.50 

161 ICSAVE ( I ) =IC ( 1 ) 

IF(NMC.EQ.l) GO TO 5 

CALL STATX(1,7,XMEAN,XC0V.XIN) 

REWIND 10 

READ (10) NOMC » IFINOMC.LT.NMO NMC=NOMC 
WRITE (6. 101 > NOMC 

101 F0RMATU7) 

READ(IO) ANOM 

WRI TE (6. 1 02 ) ANOM 

102 F ORMAT ( BE 1 6. 8) 

DO 1000 JMC=1,NMC 
DO 162 1=1.200 

162 C ( I > =CSAVE ( I ) 

DO 163 1=1,50 

163 I C ( I > =ICSAVE ( I ) 

READ(IO) AO.FRO.TBIAS.DVMOI 

IMC=JMC 

CALL 8URNA 

DVTPMA=DVTR I M 

DTAUF = 0,0 

DHPF = 0 . 0 

TB I ASF=0 ■ 0 

CALL RITE I IC 4 7) ) 

IF (IC(ll) • EQ> 0 ) GO TO 10 
CALL ERRORS ( 0 > 

I C ( 27 ) =250 
CALL BURNA 
CALL ERRORS ( 1 ) 

CALL BURN3 
CALL ERRORS (2) 

CALL BURN2 
CALL ERRORS ( 3) 

CALL 8URN1 
CALL ERRORS ( ** ) 

CALL RITE(ICUO)) 

10 CONTINUE 

XIN( 1>=DVM0I 

X I N ( 2) =DVGEOM 

X I N ( 3) =DVTRMA 

XIN(A)=DVMOI *DVTPMA 

X I N ( 5) =DTAUF 

X IN (6 ) =OHPF 

X I N ( 7) =TBI ASF 

SAVE(l»IMC)=DVMOI 

DO 200 1=2,6 

SAVE 4 1 . IMC ) =XIN 1 1 ♦ 1 ) 

200 CONTINUE 

CALL STATX(2.7 » XMEAN. XCOV ♦ X IN) 

1000 CONTINUE 

CALL STATOUT (XMEAN, XCOV. XI N. SAVE) 

GO TO 1 
5 CONTINUE 
DVMOI = 0 • 

CALL BURNA 
CALL RI TE 4 IC ( 7) ) 

IF4IC(m.EO.O)GO TO 2i> 

CALL ERRORS ( 0 ) 

CALL BURNA 
CALL ERRORS ( 1 ) 

CALL 9URN3 
CALL ERRORS (2) 

CALL BURN2 
CALL ERRORS (3) 

CALL BURNl 
CALL ERRORS (A) 

CALL RITE4IC410)) 

20 CONTINUE 
GO TO 1 
END 


Figure 1. - Listing of main program 

























SCASE 

AO 

PER 

TLAT 

TBI-AS 

DVHOI 

TWIN 

NREV 

I PHOTO 

IMIN 

JMIN 

KMIN 

LMIM 

NMC 

KEYW1 

KEYW2 

CONV 

XIPLOW 

XIPSTEP 

N5TEP 

MODE 
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Figure 3 . - Sample input. Note: Tape input overrides the values for AO, 

TBIAS, and DVMOI. 
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Figure 4. - Nominal 4-impulse transfer 
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Figure 5. - Retargeted 4 -impulse transfer. 
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Figure 6. - Retargeted 3 -impulse transfer 
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Figure 8. - Retargeted 1-impulse transfer and summary of actual end conditions. 
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